Spin-photovoltaic effect in quantum wires due to inter-subband transitions 
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We consider the current induced in a quantum wire by external electromagnetic radiation. The 
photocurrent is caused by the interplay of spin-orbit interaction (Rashba and Dresselhaus terms) 
and external in-plane magnetic field. We calculate this current using a Wigner functions approach 
taking into account radiation-induced transitions between transverse subbands. The magnitude and 
the direction of the current depend on the Dresselhaus and Rashba constants, strength of magnetic 
field, radiation frequency and intensity. The current can be controlled by changing some of these 
parameters. 
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I. INTRODUCTION 

The energy spectrum of free electrons in a perfect 
quantum wire without spin-orbit interaction consists of 
spin-degenerate subbands due to the transverse confine- 
ment in two directions. In each subband, the energy de- 
pends quadratically on the one-dimensional momentum. 
In the presence of an external radiation intersubband ex- 
citation probabilities are equal for states with opposite 
momentum. Therefore, there is no change in current as- 
sociated with external radiation. 

During the last 15 years there has been a great in- 
terest in theoretical and experimental investigations of 
photovoltaic effects and photoconductance in quantum 
wires (see Refs. [lliT3 and references therein). Mecha- 
nisms for pure spin current generation in 2D and ID sys- 
tems with spin-orbit interaction have also been discussed 
iSiiS. It is known that a photovoltaic effect is possible in 
quantum wires without inversion symmetry. For exam- 
ple, a photovoltaic effect in quantum wires with spatially 
dependent lateral confinement was predicted in Ref. U. 
In the present paper we consider a different spin-based 
mechanism for photovoltaic effect, which is very interest- 
ing in the context of the fast growing field of spintron- 
j(,g^5ji6ji7 rpj^g most important component of our scheme 
is the spin-orbit interaction. However, the spin-orbit in- 
teraction alone is not sufficient to generate a charge pho- 
tocurrent if the quantum wire is spatially homogeneous. 
Therefore, we consider a wire in an in-plane magnetic 
field which breaks the inversion symmetry. 

We can identify the following groups of inter-subband 
transitions that lead to a photovoltaic effect in quantum 
wires: 

1. Transitions between spin-splitted subbands with 
the same confinement quantum numbers. 

2. Transitions between subbands with different con- 
finement quantum numbers. 

The main difference between these two groups is that 
the first is generated by the magnetic field component 



of electromagnetic radiation, while transitions from the 
second group are due to the electric field component. In 
a recent paper^'^ a spin-photovoltaic effect in quantum 
wires due to transitions of the first type was considered. 
It was found that a special role in the photovoltaic ef- 
fect is played by transitions in which the direction of 
electron velocity changes. The importance of this veloc- 
ity inversion was outlined earlier, see e.g. Ref. 4, in 
studies of photoconductance. Here we consider intersub- 
band transitions of the second type. An important fea- 
ture in semiconductor-based quantum wires is that the 
spin-orbit interaction constants are different for subbands 
with different confinement quantum numbers. This pe- 
culiarity is essential in our scheme for the generation of 
photocurrent. 

The goal of this paper is to calculate the current 
in a quantum wire at zero bias voltage due to exter- 
nal radiation. The current as a function of radiation 
frequency is found numerically from coupled equations 
involving Wigner functions. The Wigner function for- 
malism has many advantages for investigating transport 
problemsSSiSi. Among them we mention the phase-space 
nature of Wigner functions which are similar to the classi- 
cal Boltzmann distribution functions. This feature makes 
possible to separate the incoming and outgoing com- 
ponents of the electron distribution at the boundaries 
which, in turn, facilitates the modeling of an ideal con- 
tact. The commonly used assumptions are: the distri- 
bution of electrons emitted in the quantum wire can be 
described by the equilibrium distribution function of the 
leads reservoirs, and all electron are collected by the leads 
reservoirs without reflection. In this work we extend the 
description of the transport dynamics to include inter- 
subband transitions due to electro-magnetic wave exci- 
tation. 

We show that the current is sensitive to many control 
parameters like, e.g., the spin-orbit coupling and external 
magnetic field. Therefore, the current can be used to 
determine materials parameters. The calculated current 
strength for a realistic set of parameters is of the order of 
0.1 uA and consequently can be measured using standard 
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FIG. 1: (Color online) Quantum wire with an applied mag- 
netic field in {x, y) plane, and irradiated by a electro-magnetic 
wave linearly polarized in the y-direction. 

experimental techniques. 

This paper is organized as follows. The single parti- 
cle energy spectrum and wave functions are introduced 
in Sec. m A set of coupled equations for Wigner func- 
tions is derived in Sec. IIIII The discrete model used for 
numerical solutions is presented in Sec. IIVI Numerical 
results are given and discussed in Sec. 0and concluding 
remarks are in Sec. I VII 

II. SYSTEM 



in z-direction, jiB and g* are the Bohr magneton and 
effective g-factor, Ey and w are the amplitude and fre- 
quency of electric field of the polarized electro-magnetic 
wave, and cr is the vector of the Pauli matrices. The 
effect of the external field B on the spatial motion is ne- 
glected, assuming strong confinement in the z-direction. 
The second and third terms in Q represent the Rashba 
and Dresselhaus spin-orbit interaction^^'^^ for an electron 
moving in the a;-direction, a and [3 are the corresponding 
couphng constants. 

The spin-orbit interactions included into the Hamilto- 
nian |(2Jl originate from bulk inversion asymmetry (giving 
rise to a Dresselhaus interactior^) and structural inver- 
sion asymmetry (giving rise to a Rashba interaction^i) . 
It is well known that the spin-orbit interaction con- 
stants are different for electrons in different transverse 
subbandsiSiS. In our model we assume that the Rashba 
spin-orbit interaction constant a depends on the index m 
and the Dresselhaus spin-orbit interaction constant /3 de- 
pends both on n and m, where n = 0, 1, ... and m = 0, 1, .. 
are subband indices due to confinement in y and z di- 
rections respectivelyiiSi2S In the model of rigid quantum 

wire walls P — ^ (j^im/Wz)'^ — (nm/Wy)'^^ , where 7 is a 
constant. 

At Ey ~ 0, the solutions of the Schrodinger equation 
can be written in the form 



Fig. n shows a possible experimental realization of the 
system under investigation. The two-dimensional elec- 
tron gas is split into two parts by a potential applied 
to the gate electrodes. The narrow channel between the 
gates then forms a quantum wire. Let us define a co- 
ordinate system such that the direction of the electron 
transport through the wire is in the x-direction and the 
lateral confinement is in the y-direction. We assume that 
an external magnetic field is applied in the {x, y) plane. 
Previously, several interesting investigations of quantum 
wires with spin-orbit interaction in the presence of an 
in-plane magnetic field were reportediS2iSi2iS42& 

In the quantum wire, the Hamiltonian for the conduc- 
tion electrons can be written in the form 



H = Ho + Hi 



(1) 



where 



Ho 



2m* 



-rPxCTy - -Px<7x + V{y) + U{z)-\ —a-B, 



and 



Hi 



eEyPy 



cos{ujt). 



(2) 



(3) 



Here, Hq is the time-independent part of the Hamil- 
tonian, Hi describes the interaction with the electro- 
magnetic field, p is the momentum of the electron, m* is 
the effective mass, V{y) is the lateral confinement poten- 
tial due to the gates, U{z) is the confinement potential 



(l)niiy)Vn{z), (4) 



where 
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9{..) is the step function, (j)m{y) and rjn^z) are the wave 
function of the transverse modes (due to the confinement 
potentials V{y) and U{z)). The eigenvalue problem can 
be solved to obtain 
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-a„fc H — B, 



-f3n 



(6) 



In this expression, em and E„ are the eigenvalues of de- 
coupled Scrodinger equations in y and z directions. In 
the experimental setup depicted in Fig. ^ the confine- 
ment in z-direction is stronger than the confinement in 
y-direction, thus we will consider i?i — i?o ^ ci — eo- 

In what follows we assume that in the absence of radia- 
tion the chemical potential in the wire is located between 
the ground (0,0) and first (1,0) transverse subbands so 
that only the ground subband is occupied by electrons, 
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FIG. 2: (Color online) Dispersion relations (two lowest spin- 
splitted subbands) calculated for (a) = 2T, By = 2T 
and (b) B^ = 2T, By = — 2T. These plots were obtained 
using the parameters values: m* = 0.067me, g* — —0.44, 
ao = 0.5 X 10-"eV m, /3o,o = 0.84 x 10-"eV m, /3i,o = 
0.36 X 10-"eV m, ei - eo = 3.64meV. 
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We emphasize that the direction of velocity changes in 
local extrema points of the spectrum. The radiation- 
induced transitions that we will consider conserve k. 
Since the positions of local extrema are different in the 
ground (0,0) and first (1,0) transverse subbands, transi- 
tions reversing the velocity direction are possible. The 
spectrum asymmetry results in asymmetric transition 
rates and as a result in a finite current at zero bias volt- 
age. 



III. WIGNER FUNCTIONS 
A. Interaction Hamiltonian 

Assuming a parabolic confinement in y direction we 
write the ground and the first excited wave functions of 
the corresponding transverse mode explicitly as 



<^o(2/) 



-1/4 



exp(- 



Try"- 



-1/4 



yexp(- 



„-,2 



(8) 



where y is the characteristic width of the quantum wire 
in the y direction. The energy gap between the ground 
and the first excited states can be estimated as ei — eg — 
{2fi^) / [m* y^) . Taking the form of the solution for the 
Schrodinger equation ([SJ) and (|SJ) we obtain the matrix 
form of the interaction Hamiltonian Q 



(9) 



and we focus our attention on radiation-induced transi- 
tions between these two transverse subbands. The energy 
spectrum in Eq. © is illustrated in Fig. |21 for various 
directions of B. Notice that the energy dispersion is dif- 
ferent (not simply shifted) for the ground and the first 
transverse subbands. Fig. |21 shows that the energy spec- 
trum is strongly asymmetric and significantly dependent 
on the magnetic field direction. We also note that the 
gaps between spin-splitted subbands are due to the mag- 
netic field. 



The electron velocity is defined by the slope of 

Ern,n,±{k) and is given by 



where g = {heEy) / {2m*yLd) is the coupling constant, 
is the ladder operator acting on the y component of the 
wave function: a^4'm{y) = 4'm±i{y) and s^^ is an oper- 
ator in the space of spin degree of freedom. 



„oi 



(10) 



Here s± = (1/2) (1 ± exp(iA(p„)) and A(^„ = (p„^i — 
iPn,o, and ipn,m is defined in Eq. jSjl. Neglecting high- 
oscillatory terms we can write the Hamiltonian ||^ in the 



rotating wave approximation's as 



Hi — ihg (e 



(11) 
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B. Equations for the Wigner Functions 

The Liouville-von Neumann equation for the density 
operator of the electron Pmm' ,ss' {x, x' ,t) is given by 

zhp=[H,p], (12) 

where s, s' — ± are variables associated with spin degree 
of freedom. Henceforth we omit the indices n, n' due to 
confinement in z direction since we are interested only 
in the transitions between the ground m = and m = 
1 transverse y subbands. The Wigner function can be 



obtained by integratingSi 

Wmm',ss' {R, k,t) ^ J pram'.ss' {R, Ar, t) exp{-ikAr)dAr, 

(13) 

where the density operator is written in the new spatial 
variables R = (r + r')/2 and Ar = r — r' and k is the 
electron wave vector. 

Following the standard procedure (see, e.g.j^S) we de- 
rive a set of transport equations by neglecting non-local 
correlations in the diagonal components of the Wigner 
functions 
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g [(5+1^21,+-^ + C.C.) + {S-W2l,-+ + c.c.)] , (14) 

: g [{s.W2i^+- + c.c.) + {S+W21,— + c.c.)] , (15) 
: -g [{s^W2i,^+ + c.c.) + [s+W2i^— + c.c.)] , (16) 
-g [(s+M^2i,++ + c.c.) + (s_M/2i,+- + c.c.)] . (17) 
I 

For the amplitudes of the off-diagonal components we de- fine W^^, (t) = Wmm',ss' (t) exp(— icji), and we obtain 
I 

w' vi.+ + V2,+ dW{2^++ eE, dW[^^++ 

^12.++ + ^ — + — *K2,++-w)W^i2,++ (18) 

9 [S+ {W22^++ - 1^11,++) + S_ (W^22,-+ - W^ll, + -)] , 
W"l2, + - + 7, + — T^I *(^12, + - - t^)VKi2, + _ (19) 

9 [s+ (W^22.+- - W^n.+-) + s- {W22.— - W^ii,++)] , 

W[2,- + """o^"'" + "^"—i^ - - ^)M^12,- (20) 

9 b+ {W22.— - W^ii,— ) + {W'^2,+- - W"n,-+)] : 

W^12,-+ + 7, T^^— —I. l(^12,-+ - ^)W^12,-+ (21) 

9 [s+ (W^22,-+ - W^n,-+) + (W^22,++ - W^ii,— )] , 

^11,-++ • o • S:— + ^ «(^ii,-+-^)W^ii,-+ (22) 

.9 [(s^W-2i,++ + 5+W^2i,-+) + + S-W-1'2,— )] , 

W"22,-+ + ^ Yx + *(^22,-+ - C^)VK22,_+ (23) 

= -5 [(slW^l'2,++ + <W-l'2,- + ) + (S+W^21.-+ + 5-W^21. — )] ■ 
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For completeness we also included a static electric field 
Ex along the x direction, for example, due to a bias volt- 



age. 



mm' ,ss' 



h~\E„ 



The left part of 



the equations 114I23|I describes the ballistic transport of 
the electron in the quantum wire, the right part is re- 
sponsible for excitations induced by the radiation. We 
consider the "ideal" contact boundary conditions for a 
wire of length L 



Wrr 



W, 



mm.ss 



(24) 



(25) 



where /(fc,/i) = 1/(1 + cxp[{E„„s{k) - f^)/{kBT)]) is the 
Fermi function, fc^ is the Boltzmann constant, T is the 
electron temperature and /x; are chemical potentials of 
the left/right lead, respectively. We also assume that 
only internal part of the quantum wire is irradiated, that 
is -Ey = for cc < and x > L. 

The electron charge density, electric charge and spin 
currents can be obtained from Wigner functions as 



w, 



mm.ss 



(x, k)dk^ 



(26) 



/OO 
Vm.s{k)Wmm,ss{x,k)dk, (27) 



and 



1 

i {'^^■A'^lY^^n,s)Vm.s{k) 

Wm„i^ss{x,k)dk, (28) 



respectively. Here 7 = (a;, y, z) and different ma- 
trix elements can be found in accordance with (O 



as (*„,±|cra;|*m,±) = ±COs(<Pm), (^'m,±|0-y|*m,±) = 

=F sin((/?m)- The details of the numerical model and solu- 
tion for the derived system of equations are given in the 
next sections. 



IV. DISCRETE MODEL 



The form of Eas. (|14l23|l does not allow us to solve 
the problem analytically even for the stationary case 
dW/dt = and unbiased channel Ej. = 0. The so- 
lution is complicated by different inter-subband transi- 
tions with the change of spin state s originating from 
the subband asymmetry in k domain. This effect plays 
a central role in the electric current generation and 
should be taken into account. Thus, the system of equa- 
tions H14I23|I was solved numerically for y = 25 nm, 

= 0.0001 eV and T = Q.IK. We model the domain 
X e [0, L] and k G [—k„iaxikmax] with the mesh sizes of 
Ax — L/{Nx — l) and Afc — 2kmax/{Nk — l), respectively. 
In the calculation we used = 25 and A^^ = 80 and the 
length of the quantum wire L = 2.5 fiia. The value 
kmax — 3.37 • lO^m"^ was chosen to ensure that all filled 
states in k space are taken into account. We fixed the 
values of diagonal components Wmm,ss on the boundaries 
at a; = for k with Wm,s(fc) > and at x = L for k with 
Vm,s{k) < 0, accordingly to (|24I25|I . Similarly, for off- 
diagonal components we fixed the values Wm^m' .s^s' = 
at a; = for k with (wm,s(^) + i'm'.s'(fc)) > and at x = L 
for k with {vm,s{k) +Vm',s'{k)) < 0." 

The first order upwind difference scheme was used for 
the propagation of the Wigner functions in {x, k} domain 
and second order two-step Lax-Wendroff scheme^^ to de- 
scribe the time-dependent inter-subband transfer due to 
due to interaction with the electro-magnetic field. The 
discretized Liouville equation for the Wigner function can 
be written as 



W{xi,kj,ti+i/2) = W{xi,kj,ti) 



W{x^,kj,ti+i) = W{xt,kj,ti) - At 



At r ,AW(x^,kT,ti) , 

y^k,) " ~F{W{x„k„tl)) 

v{k,)^^^^^^-F{W{x,,k,^ti+,/,)) 



r 



(29) 



where we consider the case E^ = Q and do not show ference is chosen in accordance with 
subband indices for brevity. The upwind/downwind dif- 



AW{x^,kj,ti) = < 



' W{x^,kj,ti) - W{xi^i,kj,ti) 
if v{kj) > 0, 
W{xi+i,kj,ti) - W{xi,kj,ti) 
if v{kj) < 0. 



(30) 



6 





(0,1,-) 






(0,1, +)\ / 


► 

X 



kj 




(0,0,-) 








► 




x=0 


^ x=L ^ 



FIG. 3: (Color online) The different subband transition due 
to the interaction with electro-magnetic wave are sketched. 
The interaction causes some electrons to change the direction 
of the velocity. 



The function F{W{xi, kj ,ti)) embodies the remaining 
part of Eqs. (|14I23|) which depends only on the Wigner 
function W{xi,kj,ti) itself and does not contain par- 
tial derivatives. The upwind differencing is stable'^'^ pro- 
vided the time step is small enough: Ax /At < Vmaxi 
where Vmax is the maximum possible absolute value of 
the velocity. Additionally, the time step must be much 
smaller than the highest frequency of the solution. This 
condition is satisfied by fi^^^Ai ^ 1 where r2^°^ — 

max I yj {LOmm>,ss> — w)^ -l- 1 is the maximum possible 
Rabi frequency involved in the problem. The calculation 
proved to be stable if these two conditions are met. 

The investigation of the effects of the external bias and 
charge redistribution in the quantum wire is not in the 
scope of this paper. However, the way to include these 
effects is straightforward. To consider the effect of volt- 
age applied to the quantum wire it suffices to add the x 
component of the electric field Ex{x, V{t)) in the discrete 
model (|29|l . This electric field is a function of applied po- 
tential difference V{t) = {fir — fJ'i)/^ which can be time- 
dependent and x-dependent. The latter is defined by the 
leads geometry. Additional components of electric field 
E^{n{x)) can be calculated self -consistently at each time 
step to incorporate the effect of charge redistribution. 

As was mentioned earlier, transitions between sub- 
bands can force the electron to change the direction of 
the velocity. As a result the flows of electrons moving 
in the opposite directions inside the quantum wire in- 
termix as shown in Fig. |31 Without reflections in the 
wire, the steady-state solution can be obtained by sim- 
ply imposing the condition dW/dt — and advanc- 
ing from the given values at left/right boundaries for 
v{k) ^ k. However, in the presence of radiation, the 
electron distribution at the boundary acquires an addi- 
tional component, due to the electrons whose velocity di- 



rection is changed, and the described above scheme fails. 
In order to achieve the steady-state solution we consid- 
ered the temporal dynamics of the system evolving from 
some initial state until it reached stationary conditions: 
dW/dt and I{x) ^ I{x') for aU x,x' G [0,L]. As 
a initial state we took the values of Wigner function at 
equilibrium Wmm',ss'ix,k) = Smm'Sss' f{k, fi). This cor- 
responds to the uniformly distributed electron density 
along the channel and is a solution of Eqs. H14I23|I in the 
absence of radiation and external bias: = Ey — 0. 
The chosen method of obtaining the steady-state solu- 
tion provides us also with the transient behavior and, 
thus, gives more insight into the problem. The drawback 
is a serious computational effort. The electron distribu- 
tion reaches the stationary state within the effective time 
of flight through the quantum wire. Electrons constantly 
change the direction of velocity due to the interaction 
with electro-magnetic wave, therefore this time can be 
very long compared to the time step At. Fortunately, 
the slower electrons give the smaller contribution to the 
current and a steady-state solution can be always found 
within certain accuracy. The number of steps in time 
can reach values as much as Nt ~ 10^. The results of 
the numerical simulation are presented in the following 
section [3 



V. RESULTS AND DISCUSSION 

The photo-induced current through the wire is shown 
in Fig. ^ as a function of the photon energy. The di- 
rections and strength of the magnetic field in this graph 
are the same as in Fig. 12 The amplitude of the elec- 
tric field Ey — 200 V/m used in our calculations was se- 
lected close to the electric field amplitude used in recent 
experiments34. Fig. 0] clearly shows a number of current 
peaks corresponding to different transitions. These peaks 
depend on the magnetic field strength and direction, as 
a consequence of the magnetic field dependence of the 
energy spectrum. 

In order to understand transitions leading to a specific 
peak formation, we consider in detail the current depen- 
dence on photon energy for B^- = 2T, By = — 2T (Fig. 
0Jb)). It follows from Fig. 2Jb) that the current as a 
function of photon energy has a well pronounced positive 
peak at ?ia; ~ 3.64meV, a double negative peak with a 
minimum at fiw ~ 3.67meV, a small negative peak at 
huj ~ 3.565meV, and a broad positive peak of small am- 
plitude at hu! ~ 3.75meV. Fig. represents a graphical 
determination of relevant transitions. 

The energy spectrum of two lowest spin-splitted sub- 
bands is plotted in the top panel of Fig. \E\ The bottom 
panel of Fig. O shows the energy difference between dif- 
ferent transverse subbands. By plotting horizontal lines 
corresponding to the peak energies, in the bottom panel 
and, by drawing vertical lines through the intersection 
points of those horizontal lines with energy difference, 
we can finallyidentify the points in the top panel corre- 



7 



7x10 

6x10'' 

5x10"' 



■g 4x10"" H 
^ 3x10'' 
^ 2x10'' 



1x10' - 


-1x10'" 



(a) 



'>-T-T-T-T-T-T 



3.5 



3.6 3.7 

^ft), meV 



3.8 




3.6 3.7 

hco, meV 

FIG. 4: (Color online) Current through the wire as a function 
of photon energy hu for (a) Bx = 2T, By = 2T and (b) 

Bx = 2T, iJy = 

Ey = 200 V/m. 



-2T. The parameters values are as in Fig. |5| 



sponding to the peak formation. As it was mentioned 
above, the important transitions are those that lead to a 
change of the electron velocity direction. These transi- 
tions are shown by arrows in the top panel of Fig. O 

In particular, let us consider the large negative peak 
in the current at Hu ~ 3.67meV (Fig. EJb)). Fig. 
shows that the horizontal line 3.67meV intersects only 
_Bi^o._ — £-0.0,- curve in two points. We remind that the 
electron velocity is determined by the slope of Em.n,±{k) 
according to Eq. (TJ . The left intersection point gives k- 
vector of transition in which the electron velocity direc- 
tion is conserved), because the slopes of i?i,o,- and i?o,o,- 
at this value of k are in the same direction. The right 
intersection point of 3.67me'V line with -Ei,o,- — -E'0,0,- 
curve gives a transition with a change of electron veloc- 
ity direction, specifically, with a back-scattering of left- 
moving electrons. Consequently, the electron flux from 
the right to the left decreases and, because of the negative 
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FIG. 5: (Color online) Transitions in the regions of photocur- 
rent peaks for Bx ~ 2T, By — —2T. The energy spectrum and 
difference energies are shown in the top and bottom panels, 
respectively. Horizontal lines in the bottom panel correspond 
to excitation energies in peak regions, vertical arrows in the 
top panel denote transitions when electron velocity direction 
changes. 



electron charge, a negative current appears. 

Similarly, one can consider transitions at other radi- 
ation frequencies. An interesting situation occurs for 
huj ~ 3.64meV excitation, since at this particular fre- 
quency three out of four transitions are characterized by 
the reverse of electron velocity direction. We summarize 
transitions contributing to the photocurrent at selected 
radiation frequencies in Table ^ The same analysis can 
also be applied to the result presented in Fig. ^a), but, 
because of a more distorted spectrum, the roles of dif- 



Excitation energy, meV 


Transitions 


3.565 




Eo,o,+ 


£-1,0,- 


3.64 




i?i,0,- 


; £-0,0,+ £'1,0,+ 


3.67 




Eo,o,- 


^ -Ei.o,- 


3.75 




Eo.o,- 


^ Eifi,+ 



TABLE I: Transitions giving contribution to photocurrent at 
selected radiation frequencies for Bx = 2T, By = — 2T. 
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3.6 3.7 

TiQ), meV 

FIG. 6: (Color online) Spin current components for = 2T, 
0. 



ferent possible transitions are more difRcult to interpret.. 
Moreover, we would like to note that, generally, transi- 
tion probabilities from i?o,o,- to i?i.o,+ and from -Eo,o,+ 
to Ei Q _ subbands are smaller than transition proba- 
bilities from -Eq.o.- to i?i,o.- and from -Eo.o,-!- to i?i,o,+ 
subbands because of the different spin direction in initial 
and final states. This results in a smaller current peaks 
at ?ia; ~ 3.565meV and hhj ~ 3.75meV in Fig. ^b). 

As electrons carry spin as well as charge, the external 
radiation also changes the spin current through the wire. 
Notice that even without the radiation, the spin current 
is not zero, due to the spin-orbit interaction. Fig. El 
shows spin current components at a; = for = 2T, 
By = — 2T. We note that the spin current dependence on 
the radiation frequency has features similar to the charge 
current (Fig. Elb)). However, we found that radiation- 
induced changes in spin current are much less than the 
equilibrium spin current in the wire. From an experimen- 
tal point of view, the spin currents are not so easy de- 



tectable. Therefore, the observation of this photovoltaic 
effect through spin current seems unpractical. 



VI. CONCLUSIONS 

In summary, in this paper we have investigated the 
photovoltaic effect in quantum wires with spin-orbit in- 
teraction and in-plane magnetic field. We have found 
that the peculiarities of the energy spectrum lead to a 
photocurrent generation. The dependence of the pho- 
toinduced current on the excitation frequency was calcu- 
lated numerically using the Wigner functions formalism. 
A system of coupled equations for the Wigner functions 
was derived and solved numerically for "ideal" contact 
boundary conditions. We used the first order upwind 
differencing for the propagation in the spatial domain 
and the second order two-step Lax-Wendroff differenc- 
ing for time-dependent inter-subband transitions due to 
electro-magnetic wave excitation. Stable numerical solu- 
tions were found under appropriate choices of the time- 
step Ai. The calculations can be extended to introduce 
the effects of an external bias and self-consistent poten- 
tials due to charge density redistributions which can be 
a topic for a future investigation. The frequency depen- 
dence of the photoinduced current consists of a set of 
peaks related to transitions between different points of 
the spectrum. Therefore, the energy spectrum can be re- 
constructed from photocurrent measurements. Material 
parameters, such as spin-orbit coupling constants, can be 
obtained from the analysis of the photocurrent. 
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